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I. INTRODUCTION 

In the standard model (SM) of fundamental particles, 
CP violation in weak interactions is allowed by a sin- 
gle, irreducible phase in the 3x3 Cabibbo-Kobayashi- 
Maskawa (CKM) quark flavor-mixing matrix [1, 2]. The 
unitarity of the CKM matrix, V, implies a set of rela- 
tions among its elements Vij, in particular the condition 
VJV& + V cd V c * b + V u V t t = 0, which can be depicted in 
the complex plane as a "unitarity" triangle, whose sides 
and angles are related to the magnitudes and phases of 
the six elements Vid and Vn,, where i = u,c,t. Over- 
constraining the unitarity triangle by means of precise 
measurements of all its sides and angles allows tests of 
whether the CKM mechanism is the correct description 
of CP violation. Any inconsistencies among the various 
experimental constraints would reveal effects of physics 
beyond the standard model. 

After a decade of successful operation and a total of 
about 1.3 billion BB pairs collected by the BABAR and 
Belle experiments, the three CKM angles have been mea- 
sured with varied precision. The angle (3 has been mea- 
sured with the highest precision, to around 1°, using 
B° —> (cc)K^° decays. Using a variety of two-body B 
decays (B — > irir, pir, pp and ai(1260)7r) the angle a has 
been measured to a precision of around 4° . The angle 7 
has a relatively large uncertainty, around 14°, compared 
with a and (3. The lack of precision in our knowledge of 
7 reflects the difficulty in measuring this angle. The un- 
certainties of the CKM angles quoted in this paragraph 
are taken from [3]. 

Several techniques for measuring 7 in a theoretically 
clean way are based on B meson decays to open-charm 
final states, D^°X S and D^°X S (X s = K^*, K^°). 
In these decays, the interference between the b — > cus 
and 6 —> ucs tree amplitudes, when the D° and D° de- 
cay to a common final state, leads to observables that 
depend on the relative weak phase 7. The size of the 
interference also depends on the magnitude of the ratio 
Tb and the relative phase strong phase 8b of the two 
amplitudes, which can not be precisely calculated from 
theory. They can be extracted directly from data by 
simultaneously reconstructing several related B — > DK 
decays. Many methods have been proposed to extract 
7 from B decays using D^K^* and D(*)jf(*)° final 
states (here and in the following D refers to any admix- 
ture of the neutral D° meson and its CP-conjugate D°). 
The three methods that have been used most produc- 
tively to date are the "GLW" method [4, 5], based on 
Cabibbo-supprcssed D decays to CP-eigenstates, such 
as K+R- or K^n°; the "ADS" method [6, 7], where 
the D is reconstructed in Cabibbo-favored and doubly- 
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Cabibbo-suppressed final states such as K*7r T ; and the 
"GGSZ" method [8], which studies the Dalitz-plot distri- 
bution of the products of D decays to multi-body self- 
conjugate final states, such as Ksir + ir~. A common 
problem with these methods is the small overall branch- 
ing fraction of these decays ranging from 5 x 10 -6 to 
5 x 10 -9 . Therefore a precise determination of 7 requires 
a very large data sample. BABAR has published several 7 
related measurements: GLW analyses of B* — > DK* [9], 
D*K* [10] and DK** [11] decays; ADS analyses of 
B ± _^ D (*) K ± [ 12) 13] ; DK** [11] and B° -> DK*° 
[14]; and GGSZ analyses of B* -> D^K^* [15, 16] 
and B° -> DK*° decays [17]. To date, the single most 
precise experimental determination of 7 from BABAR is 
7 = (68 ± 14 ± 4 ± 3)° and 39° < 7 < 98°, obtained from 
the GGSZ analysis of B* -> D^K^* decays [16]. In 
this measurement, the first error represents the statisti- 
cal uncertainty, the second is the experimental systematic 
uncertainty, and the third reflects the uncertainty on the 
description of the D Dalitz-plot distributions. 



II. GLW ANALYSIS OFB-} DK DECAYS 

In this paper we present the update of the GLW analy- 
sis of B* — > DK* decays based on the full BABAR dataset 
collected near the T(45) resonance. In addition to a 22% 
increase in statistics of the data sample, this study bene- 
fits from other significant improvements compared to our 
previous result [9]: 

• More refined charged track reconstruction and par- 
ticle identification algorithms, with higher purity 
and efficiency, have been employed; 

• The event shape variable J-", used to discriminate 
the signal from the continuum e + e~ — > qq back- 
ground (described in detail in Section IV) has been 
removed from the selection criteria and has in- 
stead been included in the final fit to the selected 
B candidates. This allows us to increase the sig- 
nal efficiency by about 30%. At the same time it 
provides a larger sample of continuum background 
events, thus allowing for the determination of the 
background properties directly from data (see Sec- 
tion V); 

• Better kaon/pion separation, which is needed to 
distinguish B* — > DK* candidates from the twelve 
times more abundant B — > Dir* decays, is 
achieved through the use of a global likelihood 
based not only on the Chcrcnkov angle 9c recon- 
structed by the Cherenkov detector, but also on the 
specific energy loss dE/dx measured by the track- 
ing devices. The inclusion of dE/dx in the like- 
lihood increases the kaon identification efficiency 
and decreases the pion misidentification both at 
low momentum and outside of the geometrical ac- 
ceptance of the Cherenkov detector (which is 10% 
lower than the acceptance of the tracking devices). 



In order to determine 7 from i? ± — » DK^ decays 
with the GLW method, we measure the two direct-CP- 
violating partial decay rate asymmetries, 



■A-cpk = 



_ T(B- 



D 



cp±J 



-T(B- 



D, 



CP± 



K+) 



T(B- -> D CP ±K-) ± T(B+ -> D CP± K+) 



(1) 

and the two ratios of charge averaged partial rates using 
D decays to CP and flavor eigenstates, 

R _ T(B- -» Dcp±A^-) + T(g+ -» J P C p±A'+) 
CP± ~ T(B- -> £>°A' -) + L(B+ -> fJ°A+) ' 

(2) 
where P>cp± refer to the CP eigenstates of the I? meson 
system. We then extract 7, together with the other two 
unknowns rg and 5g, by means of a frequentist proce- 
dure, which exploits the following relations [4, 5], neglect- 
ing D°-D° mixing [18]: 



Rcp± — 
Acp± = 



1 + r B ± 2tb cos 5b cos 7 , 

±2rs sin 5b sin 7 
1 ± r^ ± 2r^ cos 5b cos 7 



(3) 

(4) 



Here, r B = \A(B-^D°K-)/A(B-^D Q K-)\ is the mag- 
nitude of the ratio of the amplitudes for B~^D°K~ and 
B~—>D°K~ and 5b the difference of their strong phases. 
Taking into account the CKM factor (|T4bKs/K&K«i| ~ 
0.4) and color-suppression of the B~^-D°K~ amplitude, 
rs is expected to be around 0.1. The current world 
averages for the B ± — > DK ± GLW observables from 
the measurements in [9, 19, 20] are summarized in Ta- 
ble I. The world averages for the parameters rs and 5b 



TABLE I: World averages at 68% confidence level [21] for the 
GLW observables in B — > DK decays. 
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+ 1 
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1.10 ±0.09 0.24 ±0.07 
1.06 ±0.10 -0.10 ±0.08 



-0 104+ 0015 

'^-0.025 



and <5 B =(117ii,4) at 68% confidence 



are tb - 
level (CL) [3] 

To reduce the systematic uncertainties from branch- 
ing fractions and reconstruction efficiencies of different 
D channels appearing in the numerator and denomina- 
tor of Eq. 2, we approximate Rcp± with the double ratios 



Rcp± 
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K/tt 



R 



K/tt 



(5) 



where 
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and 



R 



Kj 71 



T(B~ -> D°K-) + T(B+ -> D°K+) 
T(B- -> D a ir-) + T(B+ -> D°tt+) 



(6) 



(7) 



Equation 5 would be exact in the limit in which the 
Cabibbo-suppressed contributions to the B ± — > Z?7r ± 
amplitudes vanish, as well as terms proportional to 
tbtd ~ 5 x 10~ 3 , as we will discuss in Section VII. This 
approximation results in a systematic uncertainty on the 
final values of Rcp± ■ 

The paper is organized as follows. In Section III we 
describe the data sample used for these measurements 
and the main features of the BABAR detector and of the 
PEP-II storage rings. In Section IV we summarize the 
procedure adopted to select B^ — > D/i ± candidates and 
suppress the main backgrounds. In Section V we intro- 
duce the simultaneous extended maximum likelihood fit 
used to extract the observables Rcp± and Acp±- In Sec- 
tion VI we explain how, by applying the same fit pro- 
cedure to selected control samples, we estimate the ir- 
reducible background present in the final samples. A 
discussion of the sources of systematic uncertainties and 
the evaluation of the uncertainties is presented in Sec- 
tion VII. Section VIII lists the final results on the GLW 
observables Rcp± and Acp±, including statistical and 
systematic uncertainties. It also contains a description 
of the statistical method used to construct frequentist 
confidence intervals for the parameters 7, 5b, and tb- 
Section IX gives a summary of our results. 



III. DATA SAMPLE AND DETECTOR 

The measurements presented in this paper use the en- 
tire BB data sample collected with the BABAR detector 
at the PEP-II asymmetric-energy B factory at the SLAC 
National Accelerator Laboratory. The BB pairs are pro- 
duced from the decays of T{AS) mesons that originate 
in collisions of 9.0 GeV electrons and 3.1 GeV positrons 
(V^=10.58GeV=M r(4S) c 2 ). In total, (467 ± 5) x 10 6 
BB pairs, approximately equally divided into B°B° and 
B + B~, have been collected in the years from 1999 un- 
til early 2008. The B meson pairs are produced almost 
at rest in the T(4S) center-of-mass (CM) frame, but the 
asymmetric beam energies boost them in the laboratory 
frame by (/3~f)cM ~ 0.56. 

The BABAR detector is described in detail else- 
where [22]. Primary and secondary vertex reconstruc- 
tion and charged-particle tracking are provided by a five 
layer double-sided silicon vertex tracker and a 40 layer 
drift chamber. Charged particle identification (PID) is 
provided by measurement of specific ionization energy 
loss in the tracking devices and of the Cherenkov radia- 
tion cone in a ring-imaging detector. Photons and elec- 
trons are identified by combining the information from 
the tracking devices and the energy deposits in the elec- 
tromagnetic calorimeter, which consists of 6580 thallium- 
doped Csl crystals. These systems are located inside a 
1.5 T solcnoidal superconducting magnet. Finally, the 
flux return of the magnet is instrumented with resistive 
plate chambers and limited streamer tubes in order to 
discriminate muons from pions. We use the GEANT4 [23] 



software toolkit to simulate interactions of particles in 

the detector, taking into account the varying accelerator 
and detector conditions. 



IV. EVENT SELECTION 

We reconstruct B^ — >• Dh^ 1 decays, where the charged 
track h is either a kaon or a pion. Neutral D mesons 
are reconstructed in the CP-even eigenstates tt~it + and 
K-K+ (D C p+), in the CP-odd eigenstates K° s n°, K° s cj) 
and KgUj (Dcp-), and in the non-CP-eigenstate K~tt + 
(D° from B~ -> D°h~) or K+n~ (D° from B+ -> 
D°h + ). CP violation in the K°-K° system is neglected, 
i.e. the K° is assumed to be a pure CP = +1 cigenstatc. 
The Dcp daughters are reconstructed in the decay modes 
K° — > 7T+7r~, <p — > K + K~ and lu -> 7r~7r + 7r . 

We optimize all our event selection requirements by 
maximizing the significance of the expected B ± — > DK^ 
signal yield, defined as N s i g /y/N s [ g + A^kg, where N s [ g 
(-^Vbkg) is the expected signal (background) yield. The 
optimization is done for each D decay channel using sim- 
ulated signal and background events, which are generated 
with the EVTGEN software package [24] . 

Neutral pions are reconstructed by combining pairs 
of photon candidates with energy deposits larger than 
30MeV that arc not matched to charged tracks and 
whose energy deposition profile is consistent with that 
expected from a photon. The photon pair invariant mass 
is required to differ from the nominal 7r° mass [25] by less 
than 2.5 times its resolution (a fa 6 MeV/c 2 ) and the total 
7T° energy in the laboratory frame must be greater than 
240 MeV for D -> K° s -k° and 210 MeV for w -> tt+tt-tt . 

Neutral kaons are reconstructed from pairs of oppo- 
sitely charged tracks with invariant mass within 2.5er 
(er w 2.1 MeV/c 2 ) of the nominal K® mass [25]. The ratio 
between the K® signed 3-dimcnsional flight length and its 
uncertainty, determined from the position of the K® and 
the D decay vertices and the K® momentum direction, 
must be greater than 1.9, 2.0, and 2.2 for D -> K°tt , 
D — > K®<f>, and D — > KgU), respectively. 

The </> candidates are reconstructed from pairs of oppo- 
sitely charged tracks passing kaon identification criteria 
with typical kaon selection efficiency of sa 98% and pion 
misidentification of w 15%. The two tracks are assigned 
the kaon mass hypothesis and their invariant mass is re- 
quired to be within 6.5 MeV/c 2 of the nominal <f> mass [25] 
(the resolution is a = 1.0 MeV/c 2 and the natural width is 
T0 = 4.3 MeV). We also require that the helicity angle Op 
between the flight direction of one of the two kaons and 
the D flight direction, in the <f> rest frame, satisfies the 
condition |cos#y| > 0.4. This requirement exploits the 
fact that in D — > K®<fi decays the <fi is produced in a lon- 
gitudinally polarized state, thus cos 0h follows a cos 2 Oh 



of the nominal u mass [25] (the resolution is a = 
6.9MeV/c 2 ). The charged pion candidates are required 
to pass pion identification criteria with pion selection 
efficiency around 98% and kaon misidentification rate 
around 12%. To improve the uj momentum resolution, 
the invariant mass of the two photons forming the ir° can- 
didate is constrained to the nominal ir° mass. We define 
On as the angle between the normal to the w decay plane 
and the D momentum in the u> rest frame, and O^n as 
the angle between the flight direction of one of the three 
pions in the u rest frame and the flight direction of one of 
the other two pions in the two-pion rest frame. The quan- 
tities cos On and cos V ^ follow cos 2 On and (1 — cos 2 V7r ) 
distributions, respectively, for the signal, and are almost 
uniformly distributed for wrongly reconstructed u) candi- 
dates. We require the product cos 2 On sin 2 T7r > 0.046. 
Neutral D candidates are formed from two-body com- 
binations of K^, n , K®, 7T°, (j) and w candidates consis- 
tent with one of the six D decay channels under study. 
To improve the Dcp- momentum resolution, the invari- 
ant masses of the 7r° and K® daughters are constrained 
to the nominal ir° and K® masses. To suppress poorly 
reconstructed D candidates and candidates from ran- 
dom combinations, we perform a geometric fit of the D 
daughters to a common origin, and reject D candidates 
for which the \ 2 probability of the vertex fit is lower 
than 0.01%. The invariant mass of a D candidate Mp 
must be within a range that corresponds to slightly more 
than twice the Mp resolution, which varies from about 
6 MeV/c 2 for the K%<$> channel to about 44 MeV/c 2 for the 
K®tt° channel. We apply the following particle identifica- 
tion criteria to the charged daughters of the D meson: in 
D — > 7r + 7r~ , the two pion candidates must pass the same 
pion identification criteria adopted in the reconstruction 
of w -> ir~Tv + ir°; in D ^ K + K~, the two kaon can- 
didates are required to pass tighter kaon identification 
criteria than those applied to the <f> daughters (typical 
kaon selection efficiency around 94%, and pion misiden- 
tification rate around 6%); in D —t K~ir + , the kaon 
candidate must pass the same kaon identification criteria 
required for the </> daughters. In order to reduce the large 
combinatorial background from random combinations of 
tracks and photons in e + e~ — > qq events (q = u, d, s, c), 
we put requirements on the cosine of the D decay angle, 
I cos Op I . We define Op as the angle between one of the D 
daughters in the D rest frame, and the direction of the D 
meson in the B rest frame. Due to angular momentum 
conservation we expect the distribution of cos#£> to be 



uniform for B~ 



D/i±,Z) 



and D -> K Q s n° sig- 



distribution, while in <j> candidates not from D — > K®(j) 
decays, cos Op is approximately uniformly distributed. 

The ui candidates are reconstructed from ir + ir~~ir° com- 
binations with invariant mass within 17 MeV/ c 2 (2r w ) 



nal events, while for qq events the distribution is strongly 
peaked at ±1. We require | cos#£i| < 0.74 (0.99) for the 
B ± -4 Dh ± , D -> 7T+7T- (D -> #°7T°) channel. 

The invariant mass distributions of the reconstructed 
D candidates, after all the other selection criteria de- 
scribed in this section have been applied, are shown in 
Fig. 1. 

We reconstruct B meson candidates by combining a 
neutral D candidate with a track h^. For the D^-Kir 
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FIG. 1: Distributions of the difference between the D candidate's invariant mass and the nominal D mass [25], as measured 
in the B — ► Dh samples. All selection criteria described in Section IV, except that on the D invariant mass Mo, have been 
applied, including the x 2_ based candidate selection. In addition we reduce the background by requiring the fit variables to 
satisfy tties > 5.27GeV/c 2 , AE > —0.05 GeV, and T > —0.25. The AMo selection requirements are depicted by the vertical 
lines. 



mode, the charge of the track h must match that of the 
kaon from the D meson decay. This selects b — > c medi- 
ated B decays B~ ->■ D°h,- and B+ -> D°h+. The con- 
tamination from b — > u mediated B decays followed by 
doubly-Cabibbo-suppresscd D decay, i.e. B~ — > D°K~, 
D° — > K~tt + , and from D°-D° mixing is negligible. 
In the S ± — > Dh^, D — > tt + tt~ channel we require 



tion, the neutral D invariant mass is constrained to the 
nominal D° mass [25] for all D decay channels. 

We identify signal B — > DK and B —* Dtt candidates 
using two kinematic variables: the difference between the 
CM energy of the B meson (Eg) and the beam energy, 



AE = E* B - V^/2 , 

that the invariant mass of the (/i±7tF) syste mjs greater and thc beam-energy-substituted mass, 
than 1.9 GeV/ <r to reject background from B —tDw , 



D°^K- 



and B-^K 



K* U ->K 7T+ decays and 



their CP conjugates. Here 7r is the pion from the D and 
h is the track from the B candidate taken with the kaon 
mass hypothesis. To improve the B momentum resolu- 



"T-ES = J (s/2 + Pee • p B ) 2 /E% 



■Pb, 



(8) 



(9) 



where (Eb,Pb) and (E ee ,p ee ) are the four-momenta of 
the B meson and of the initial e + e~ system, respectively, 
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measured in the laboratory frame. The toes distributions 
for B^ — > Dh^ signals are centered at the B mass [25], 
have a root-mean-square of approximately 2.6MeV/c 2 , 
and do not depend strongly on either the D decay mode 
or the nature of the track h. In contrast, the AE dis- 
tributions depend on the mass assigned to the track h. 
We evaluate AE with the kaon mass hypothesis so that 
the peaks of the distributions are centered near zero for 
B ± — » DK^ events and are shifted by approximately 



+50MeV for B d 



Z?7r ± events. The AE resolution 



depends on the kinematics of the decay, and is typically 
16MeV for all D decay modes under study after the D 
invariant mass is constrained to its nominal value. We re- 
tain B candidates with toes and AE within the intervals 
5.20 < toes < 5.29GeV/c 2 and -80 < AE < 120 MeV, 
which define the region for the fit described later. 

In order to discriminate the signal from e + e~ — > qq 
background events, denoted qq in the following, we con- 
struct a Fisher discriminant J- based on the four event- 
shape quantities Lf OE , |cos#^|, \cosB* B \ and -ff^ OE . 
These quantities, evaluated in the CM frame, are defined 
as: 

• ^20° E = L2/L0 is the ratio of the second and ze- 
roth event shape moments of the energy flow in the 
rest of event (ROE), i.e. considering all the charged 
tracks and neutral clusters in the event that are not 
used to reconstruct the B candidate. They are de- 
fined as L 2 = J2i Pi cos2 ®i an d L = J2i Pi > where 
Pi are the momenta and 0i the angles of the charged 
and neutral particles in the ROE, with respect to 
the thrust axis of the B candidate's decay products. 
The thrust axis is defined as the direction that max- 
imizes the sum of the longitudinal momenta of the 
particles used to define it; 

• 9?p is the angle between the thrust axis of the B 
candidate's decay products and the beam axis; 

• 9* B is the angle between the B candidate momen- 
tum and the beam axis; 

• ^2o° E = H2/H0 is the ratio of the second and ze- 
roth Fox- Wolfram moments H2 and Hq [26], com- 
puted using charged tracks and photons in the 
ROE. 

The quantity J 7 is a linear combination of the four afore- 
mentioned event-shape variables: 



For events with multiple B ± —$■ Dh^ candidates 
(about 16% of the selected events), we choose the B can- 
didate with the smallest x 2 = E c ( M c - A'/ c PDG ) 2 /(er^ /c + 
r 2 ) formed from the measured and true masses, M c and 
A1^ BG , of all the unstable particles c produced in the 
B decay tree (D, ir°, Kg, </>, u>), scaled by the sum in 
quadrature of the resolution UM a of the reconstructed 
mass and the intrinsic width T c . From simulated signal 
events, we find that this algorithm has a probability to 
select the correct candidate between 98.2% and 99.9% 
depending on the D decay mode. We also find that the 
algorithm has negligible effect on the Mo distributions. 

We compare the distribution of each selection variable 
in data and simulated events after the requirements on 
all other variables have been applied. In order not to 
introduce biases that may artificially enhance the signal 
yield, we perform a blind study by explicitly removing, in 
this comparison, events consistent with the B ± — > DK^ 
signal, i.e. those with |toes — TOb| < 10 MeV/c 2 , \AE\ < 
40 MeV, T > —0.8 and track h passing kaon identi- 
fication criteria. We find excellent agreement between 
data and simulated events, both for events consistent 
with the B ± -*■ Dtt* signal (|to es - m B \ < 10 MeV/c 2 , 
\AE - 50 MeV| < 40 MeV, T > -0.8 and track h failing 
the kaon identification criteria) and for background-like 
events. We correct for small differences in the means and 
widths of the distributions of the invariant masses of the 
unstable particles and of toes and AE both when apply- 
ing to data the selection criteria obtained from simulated 
events and in the final fit described in the next section. 

The total reconstruction efficiencies, based on simu- 
lated B ± — > DK^ events, are summarized in the sec- 
ond column of Table II. For the reasons explained in 



TABLE II: Reconstruction efficiency for B — > DK from sim- 
ulated events. We also quote the efficiency and purity in a 
signal-enriched subsample (see text for details). 



D° mode 


Efficiency after 


Efficiency in 


Purity in 




full selection 


signal-enriched 


signal-enriched 






subsample 


subsample 


K-TT + 


52% 


22% 


96% 


K+K~ 


44% 


18% 


85% 


TT + TV~ 


38% 


17% 


68% 


K S TT 


24% 


10% 


83% 


kU 


20% 


9% 


91% 


K° s uj 


10% 


4% 


71% 



T = ciL^ )E + ca| cos6* T \ + c 3 \cose* B \ 



C4^ 2 R o OE . (10) 



The values of the coefficients Ci are the ones which max- 
imize the separation between simulated signal events 
and a continuum background sample provided by off- 
resonance data, taken « 40 MeV below the T(45) res- 
onance. The maximum likelihood fit described in Sec- 
tion V is restricted to events with T within the interval 
— 1.5 < J- < 1.5, to remove poorly reconstructed candi- 
dates. 



Section II, the efficiencies are 40% to 60% higher than 
in our previous study of the same decay channels [9]. 
The efficiencies obtained for B ± — s- £>7r ± events from 
the simulation are statistically consistent with those for 
B ± — > DK^, where the D meson is reconstructed in 
the same final state. For illustration purposes we define 
a signal-enriched sample for each D decay mode, con- 
taining all B^ —> Dh^ candidates satisfying the criteria 
-40 < AE < 100 MeV, 0.2 < T < 1.5, 5.275 < to es < 
5.285 GeV/c 2 , and whose daughter track h passes charged 
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kaon identification criteria. The typical kaon efficiency is 
pa 77% and the pion misidentification rate is pa 2%. The 
reconstruction efficiencies and the expected purities for 
the signal-enriched subsamples, determined on simulated 
data, are listed in Table II. 



where the iVJ are the expected yields in each component 
j. In case of negligible correlations among the fit vari- 
ables, each probability density function (PDF) V factor- 



ies as: 



V(m ES , AE, T) = V{m ES )V(AE)V(T). 



(13) 



V. MAXIMUM LIKELIHOOD FIT 

We measure R K L and Acp± using simultaneous ex- 
tended and unbinned maximum likelihood fits to the dis- 
tributions of the three variables AE, toes, an d T of B 
candidates selected in data. The dataset is split into 
24 subgroups by means of three discrete variables: the 
charge 77 = ±1 of the reconstructed B meson (x2 sub- 
groups); the two-body D decay final state X (x6), allow- 
ing for a more accurate description of the corresponding 
probability density functions compared to the larger CP± 
subgroups; and a PID variable denoting whether or not 
the track h from the B passes (p) or fails (/) charged 
kaon identification criteria (x2). The pion misidentifi- 
cation rate of these criteria is determined directly from 
data as described later, and is expected from simulation 
to be around 2%. The corresponding kaon identification 
efficiency is (77 ± 1)%, as determined from the signal 
MC samples after weighting the bidimcnsional distribu- 
tion of the momentum and polar angle of the track h by 
the ratio of the analogous distributions observed in MC 
and data kaon control samples. The uncertainty on the 
kaon identification efficiency is dominated by the system- 
atic contribution from the uncertainties on the weights. 
We perform in total three simultaneous fits to these 24 
subgroups: one fit for the two CP-even D final states 
(8 subgroups), one for the three CP-odd D final states 
(12 subgroups), and one for the D — ► Ktt decay (4 sub- 
groups). 

The likelihood function C for each of these simultane- 
ous fits has the form 



,-JV 



m 



N 1 ' 



N s 



JVsim^uAEi,^;^, (11) 



where s ranges over the subgroups under consideration, 
N s is the number of events in subgroup s, n is the total 
number of events in the fit n = ^2 S N S , and N is the 
expected number of events. We minimize — ln£ with 
respect to the set of fit parameters v specified later. 
The probability V 8 ,i = Ps(toesj, AEi,^) for an event 
i is the sum of six signal and background components: 
B ± — > DK^ signal, B ± — > Dn^ signal, background can- 
didates from e + e~ — > qq events, irreducible background 
arising from charmless B ± — > XK^ and B^ — > Xtt^ de- 
cays, and background candidates from other BB events 
(reducible BB background): 



N M V„ 



s ' s,i 



+ N?*Pff 



Nf n V 



Xtt 
sA 



N XK V 



,XK 
s.i ; 



(12) 



The irreducible BB background originates from events 
where a B meson decays to the same final state Xh as 
the signal, but without the production of an intermediate 
charmed meson in the decay chain. When exploiting the 
AE, toes, and T variables, this background is therefore 
indistinguishable from the signal. As an example, the 
decay B ± -¥ K + K~K ± (X = K + K~) is an irreducible 



background for B ± — > Dcp+K ± , Dcp+ 



K + K~ 



As described later in Section VI, the irreducible back- 
ground yield can be estimated by studying sideband re- 
gions of the D candidate invariant mass distribution, 
and can then be fixed in the final fit, where we assume 

pDh _ -pXh 

We express the signal yield parameters N® K and 
N^ through the CP asymmetries A^ K and A^ of 
B ± -> DI^, D^X and B ± -> L> 7 r ± , D^X, their 
branching fraction ratios, R^/^, the total number N®£ x 

of -B — > Dtt^ 1 , D — >• X signal events, the true kaon iden- 
tification efficiency e of the PID selector, and the pion 
misidentification rate to of the PID selector: 



N 



V, 



DK 
V,P,X 

DK 

■n,f,x 



1 



= ^(1-A) N& tX R$ /ir E 



(14) 



= 5 (1 " VAI K ) N& tX Rl h (1 - e) , (15) 



N 



Dtc 

n,v,x 



]\tDtv 



Da 



1-vAZJ N»l x m, 



1 - f^£J A& (1 - TO 



(16) 

(17) 



Because the ratios R^., are small, the fit is not able 

K/7T ' 

to determine the value of e. Therefore we fix it to the 
aforementioned value of e = (77 ± 1)%. The reconstruc- 
tion and selection efficiencies for true B ± —± DK^ and 
B ± — > Z?7r ± candidates, where the D meson decays to the 
same final state, are assumed to be identical. A system- 
atic uncertainty is assigned due to this assumption (see 
Section VII). The simultaneous fit to the two CP-even 
modes constrains 



A DK 

TfTT + TT~ 

Rr/tt 



A K+K ~ — A 
A DK = A CP+ 

tdK+R- _ p+ 
^K/tt — n K/-K 



(18) 
(19) 



while the simultaneous fit to the three CP-odd modes 
constrains 



A DK 






A DK 



A 



,K" 



DK 

K°LJ 



A C p- 



™K/n — "-K/-K — "-K/-K 



= R 



K/-K 



(20) 

(21) 



The to-es distributions of the signal components are pa- 
rameterized using an asymmetric Gaussian shape, i.e. a 
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Gaussian with different widths on both sides of the peak. 
We use the same shape for B^ — > DK^ and B ± — > Dir^, 
so the toes B —>■ DK signal shape (whose parame- 
ters are floating in the fit) will mostly be determined by 
the much more abundant B — > Dir control sample. 
Since the selection efficiencies for the two channels are 
the same, we expect the number of reconstructed candi- 
dates from B ± — > Dir^ 1 to be about twelve times higher 
than for B ± — > DK^. We have checked that the toes 
shapes for i? ± — » DK^ and B ± — > Dir^ are consistent, 
and that the assumption that they are identical does not 
bias the parameters of interest. 

The AE distribution of the B^ — > DK^ signal compo- 
nent is parameterized with a double Gaussian shape. The 
core Gaussian has a mean close to zero, a width around 
16MeV and, according to the simulation, accounts for 
about 90% of the true B ± — > DK^ candidates. The sec- 
ond Gaussian accounts for the remaining 10% of candi- 
dates whose energy has been poorly measured. The mean 
and width of the core Gaussian are directly determined 
from data, while the remaining three parameters (the dif- 
ference between the two means, the ratio between the two 
widths and the ratio of the integrals of the two Gaussian 
functions) are fixed from the simulation. In contrast to 
the toes case, the B ± — > Div^ 1 AE shape is not the same 
as for B ± — > DK ± . This is due to the fact that we always 
assign the kaon mass hypothesis to the track: the wrong 
mass assignment, in the case of S ± — > Dir^ 1 , introduces 
a shift to the reconstructed energy of the pion and thus 
to AE, since AE = E* B - y/H/2 = E* D + E* h - y/1/2. The 
shift depends on the magnitude of the momentum p of 
the track h in the laboratory frame, 



AE shitt (p) = 7CM 



m 2 K + p 



-V? 



(22) 



Therefore we parameterize the B — > Dir AE signal 
component with the sum of two Gaussians whose means 
are computed event-per-event by adding Ai? S hift (p) to 
the means of the Gaussian functions used to describe the 
B ± — > DK^ AE signal. The other parameters of the 
B ± -> Dn ± and B ± ->■ DR r± AE distributions (the two 
widths and the ratio of the integrals) are identical. Again, 
we exploit the B — > Dir control sample to determine 
the shape of the B^ — > DK^ signal. In the case of the 
high statistics flavor mode D — ¥ Kir, we add a linear 
background component to the double Gaussian shape to 
account for misreconstructed events, which peak in Toes 
but not in AE. The ratio between the integral of the 
linear component and that of the two Gaussian functions 
is fixed from simulated signal events. 

For the reducible BB background, Eq. 13 does not hold 
because of significant correlations between the AE and 
Toes distributions. This reflects the fact that this back- 
ground is composed of two categories of B candidates 
with different Toes and AE distribution: 

• B candidates formed from random combinations of 
charged tracks and neutral objects in the event, 
which populate the whole rriES-AE plane; 



B candidates from B ± -> Dp ± , B ± ->■ DK^, 
B ± -> D*h ± (D* — > Dir), where a pion from the 
p, K* or D* decay is not reconstructed. These 
candidates peak in toes close to the B mass, but 
with broader resolution compared to the signal, and 
are shifted towards negative AE values, typically 



peaking at AE 



-to^c 



2 , therefore outside of the 



AE fit region; however, the tail on the positive side 
of the distribution extends into the AE fit region. 

We parametrize the m^s-AE distribution of the BB 
background by means of two factorizing components: 

V BE (m E s,AE) = f x g pcak (mEs)h pca k(AE) + 

(1-/) x g C o„t(TOEs)/w(A£). (23) 

The toes component of the peaking part, 5 P eak("ZEs), 
is parameterized with a Gaussian function for X = 
iT+iT-, K° s uj, K°cf>. For X = K+K~ , A'°7r° we use the 
"Crystal Ball" lineshapc [27], an empirical smooth func- 
tion that better describes the non-Gaussian tail on the 
negative side of the distribution, 



C(x) 



i r^e- 1 ^ (A-\a\-i 
I exp(-is 2 ) 



— n 
X ) X < 



\a\ 
In I 



x> 

'(24) 

with x = (x — fi)/ 'a and x — > ~x for a < 0. For X = Kir 
we use an empirical function of the form: 



1 



N(x) = exp I -— g {In [1 + At(x -/jl)]+t 



2t 2 



*>)■ 



(25) 



with A — sinh(r \/ln 4) / (err Vhi 4) . Here /i is the posi- 
tion of the peak, while a and r are parameters related 
to the width of the distribution on the two sides of the 
peak. The AE component of the peaking part h pca k(AE) 
is described with a simple exponential function for the 
five CP self-conjugate D final states, and with a Landau 
function for the non-CP-eigenstate final state. The BB 
purely combinatorial background component is described 
by the 2-dimcnsional product of a linear background, 
h CO nt(AE), and an empirical function introduced by the 
ARGUS collaboration [28], g CO nt(™Es) = A(niEs/ r no)'- 



A{x) = x(l - x 2 ) p exp (-C [1 - x 2 } ) 



(26) 



where too = y / i/(2c 2 ) = 5.29GeV/c 2 is the kinematic 
endpoint of the toes distribution. All the parameters of 
the BB background toes-AE 1 distribution are fixed from 
simulated BB events. The only exception is the width of 
the Landau function used for h vea ^{AE) in X — K~ir + . 
This parameter controls the behaviour at low AE values, 
AE w —80 MeV, where we find the simulation not to be 
sufficiently precise given the high statistics of this chan- 
nel. We note that the shape parameters differ across the 
six final states, but are similar across the charge and PID 
selector subgroups belonging to one final state. 

In qq events, B candidates arise from random combina- 
tions of charged tracks and neutral particles produced in 
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the hadronization of the light quark-antiquark pairs pro- 
duced in e + e~ collisions. Similarly to the combinatorial 
component of the BB background, the qq background 
distribution in the m^s-AE plane is parameterized by 
the product of an ARGUS function in toes and a linear 
background in AE. We float the slope of the linear com- 
ponents, while the parameters of the ARGUS function 
are fixed, in each D final state, from simulated qq events. 
They are in good agreement across the final states and 
other subgroups. 

The J- distributions are parameterized in a similar way 
for all fit components. We find that the distributions of 
S ± — > DK ± and S ± — > Dit^ 1 signal events are consis- 
tent with each other, as expected since their kinematics 
are very similar, and choose to parameterize them with 
the same shape. For this we use the sum of two asym- 
metric Gaussian functions. Some channels with lower 
statistics don't require the full complexity of this param- 
eterization: in those cases we use a single asymmetric 
Gaussian, a double Gaussian, or a single Gaussian. In 
particular we use: for the signal components a double 
asymmetric Gaussian, except for X = Kg<j), where a 
double Gaussian function is adopted; for the BB back- 
ground components a double asymmetric Gaussian in 
case of X = K®ir° ,K~ir + , a double Gaussian in case 
of X = KgLo, and a single Gaussian otherwise; for the qq 
background components a double asymmetric Gaussian, 
except for X = K®(f), where we use a single Gaussian. 

In summary, the floating parameters of the fits are: 
all parameters related to the signal yields, and therefore 
to the GLW parameters, as given in Eqns. 14-17, except 
e; all background yields and CP-asymmetries except the 
irreducible background yields and asymmetries, the BB 
asymmetries for CP— modes and for the (CP+,p) sub- 
groups (B — > Dcp+h candidates where the track h passes 
the kaon identification criteria) , and the BB yield in the 
(K®(fi,p) subgroup; selected shape parameters, namely 
the overall width and mean of the AE signal, the Toes 
signal shape, and the AE and JF shape for qq background. 
A full list of the floating parameters can be found in Ta- 
ble VIII. The non-floating parameters are fixed to their 
expectations obtained from simulation or, in case of the 
irreducible background yields, to values obtained from 
data control samples (see next section). Non-floating CP 
asymmetries are fixed to zero. We assign systematic un- 
certainties due to the fixed parameters. 

We check that the fitter is correctly implemented by 
generating and fitting a large number of test datasets us- 
ing the final PDFs. In this study, we include an analytic 
description for the conditional variable Ai^ sh ift- The 
residuals for a given parameter, divided by the measured 
parameter error, should follow a Gaussian distribution 
with zero mean (p) and unitary width (a). We observe no 
significant deviations from the expected distribution. In 
particular, R K / V shows the largest shift from zero mean 
(fi = —0.06 ±0.07) and Acp+ shows the largest deviation 
from unity width (a — 1.13±0.06) among the parameters 
of interest. 



We investigate fit biases, arising from possible discrep- 
ancies between the true signal distribution and the cho- 
sen fit model, by fitting a large number of test datasets, 
in which the B ± — > DK^ and S ± — > Dir^ 1 signal com- 
ponents are taken from simulated samples of sufficient 
statistics, while the background components are ran- 
domly generated according to their PDFs. Of all floating 
parameters, only Rk/tt acquires a significant bias, result- 
ing in corrections of 0.5 and 1.0 times the expected sta- 
tistical uncertainties on these parameters in the CP and 
flavor modes, respectively. This bias is caused by small 
differences in the AE distributions of the signal compo- 
nents across the kaon PID subgroups (p and /), which 
the final PDF does not account for. A second, smaller 
contribution to this bias is a small discrepancy between 
the AE(tt) signal shape of B ± —> Dit^ 1 events and the 
AE(K) shape of B ± —> DK^ events. The biases in 
the Rk/tt parameters are correlated, and partly cancel 
in the ratio, resulting in a smaller bias on the GLW pa- 
rameters Rcp±- The largest (smallest) remaining bias is 
0.12 (0.05) times the expected statistical uncertainty for 
Rcp+ {Acp-)- We correct the final values of the param- 
eters Acp and Rk/tt for the observed biases, and assign 
systematic uncertainties to these corrections. 



VI. IRREDUCIBLE BACKGROUND 
DETERMINATION 

As discussed in the previous section, the irreducible 
background arises from charmless B ± — > Xh decays, 
which have the same final states as the B ± — > D(— > 
X)h signal and therefore the same distribution of the 
three fit variables AE, jtles, and J- ' . 

In the D° — > K~iv + flavor mode, the irreducible back- 
ground - taking into account the measured branching 
fractions for B ± ->■ K^K^ and B ± -► K^^ir* [25] 
and a selection efficiency of ss 1%, estimated from sim- 
ulated events - is negligible compared to the expected 
signal yields (about 3400 B ± -5> £>A ± and 45000 B ± -> 
Dtt expected signal events). On the other hand, in the 
CP modes, where the signal yields are expected to be an 
order of magnitude lower than in K~tt + , and the upper 
limits for the branching ratios of S ± — > Xh^ 1 decays are 
at the 10 -5 level, we cannot a priori exclude a relevant 
irreducible background contribution. 

We estimate the irreducible background yields in our 
sample by exploiting the fact that the D invariant mass 
distribution for this background is approximately uni- 
form, while for the signal it is peaked around the nom- 
inal D mass. Therefore we can select a control sample 
containing irreducible background candidates, but with 
the signal strongly suppressed, by applying the same se- 
lection as for the signal, with the only difference that the 
D invariant mass is required to lie in a region (D invari- 
ant mass sidebands) which is separated by at least a few 
(Tm d from the nominal D mass (see Table III). We then 
perform an extended maximum likelihood fit to the toes, 
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AE, and JF distributions of the control sample in order 
to measure the irreducible background yields in the D in- 
variant mass sidebands. The fit is similar to the nominal 
one described in the previous section. However, due to 
the limited statistics available in the sidebands, we are 
forced to fix more parameters compared to the nominal 
fit; in particular, we fix any possible charge asymmetry of 
the B^ — > Xh^ decays to zero (a systematic uncertainty 
is assigned to this assumption). Finally, since the D can- 
didate invariant mass distribution of the irreducible back- 
ground is approximately uniform, we scale the obtained 
yields by the ratio of the widths of the D signal and con- 
trol sideband mass regions to obtain the irreducible back- 
ground yield N xh (scale factor in Table III). Table IV 
shows the scaled irreducible background yields that enter 
the final fit. 



TABLE III: D mass sideband definitions, the scale factor de- 
fined as the ratio of the widths of the D mass signal and 
sideband regions. 



D decay 




Mb sideband regior 


i 


Scale 


mode 




(MeV/c 2 ) 




factor 


K+K- 


[1794.5 


-1834.5], [1884.5- 


1914.5] 


0.43 


n n~ 


[1814.5 


- 1839.5], [1889.5 - 


1934.5] 


0.48 


1^0 
K S 1T 


[1774.5 


-1804.5], [1924.5- 


1954.5] 


1.67 


K° S 0J 


[1794.5 


-1829.5], [1899.5- 


1934.5] 


0.69 


K° s <f> 


[1794.5 


-1834.5], [1894.5- 


1934.5] 


0.28 



TABLE IV: Irreducible background yields estimated from Md 
sidebands in data. 



D decay mode N 



N J 



K+K~ 

A S 7T 

K° s uj 



93 ± 10 


-5 ± 8 


4 ± 6 


± 9 


-4 ± 9 


65 ± 23 


3 ± 6 


± 8 


0.5 ± 0.7 


1.4 ± 1.0 



VII. SYSTEMATIC UNCERTAINTIES 

We consider nine sources of systematic uncertainty 
that may affect the GLW parameters Acp± and Rcp±- 
Their contributions are summarized in Table V. 

First, we estimate the influence of fixed parameters of 
the nominal PDF. We perform a large number of test 
fits to the data, similar to the nominal fit. In each of 
these test fits the fixed parameters are varied according 
to their covariance matrices. From the resulting distri- 
butions we calculate the systematic covariances of the fit 
parameters Acp± and Rk/k- The parameters responsi- 
ble for the largest uncertainty are the toes endpoint mo, 
and parameters related to the measured yields, e.g. bb 
background asymmetries and the efficiency of the kaon 
selector. 



The uncertainties in the irreducible background event 
yields introduce a systematic uncertainty in the B ± — > 
Dcph ± yields and therefore in Rcp±- Likewise, any 
charge asymmetry in this background would affect the 
measured values of Aqp±. We again perform a series of 
test fits to on-pcak data, where we vary the B^ — > Xh^ 
yields and asymmetries by their uncertainties. For the 
latter, we take the uncertainties to be ±10% for X = 
K + K~ and ±20% for the other modes, which are conser- 
vative estimates consistent with the existing upper limits 
on the CP asymmetries in those decays [21]. For Acp+, 
the possible CP asymmetries in the peaking background 
dominate the systematic error. 

As explained in Section V, we correct the fit results 
for biases observed in Monte Carlo studies. We take the 
associated systematic uncertainties to be half the size 
of the bias corrections, summed in quadrature with the 
statistical uncertainties on the biases. The latter are due 
to the limited number of test fits used to estimate the 
corrections. 

We investigate a potential charge asymmetry of the 
BABAR detector, due to a possible charge bias in tracking 
efficiency (e.g. K + vs K~) and/or particle identification. 
Our analysis includes a number of control samples, in 
which the CP asymmetry is expected to be negligible: 
the six B ± — > Dir^ 1 samples and the B ± — > DK^ flavor 
mode (D — > Kir). The weighted average of the charge 
asymmetry in the control samples is (—0.95 ± 0.44)%, 
from which we assign uncertainties of 1.4% to both Acp+ 
and Acp~. We consider these uncertainties to be 100% 
correlated. 

The measured CP asymmetry in B ± — > DK ± , D —> 
K®4>, can be diluted by the presence of B ± — > DK ± 
decays followed by D decays to the same final state 
K®K + K~ as the signal but with opposite CP content, 
such as D — > K®ao, a® — > K + K~ . The same can hap- 
pen in the B ± — > DK ± , D —> K®lu analysis with back- 
grounds from B ± -)■ DK ± , D -> A'°7r~7r + 7r°. This back- 
ground can also affect the ratios Rcp-- It is possible to 
obtain correction factors to both Aqp- and Rcp- from 
a fit to the distributions of the relevant helicity angles, 
cos On and cos 9r for K°oj and K°cj), respectively. The fit 
is performed on dedicated B ± — > Dir^ samples, in which 
the selection requirements on the helicity angles have not 
been applied. It can be shown [29] that for these two final 
states the observed charge asymmetries and ratios should 
be corrected by a factor 



/(true 
A CP 



ntruc 
n K/ir 



4obs 
A CP 



l±/eM 2 ^ 

i-/,M 2 ' 
bs 1 + /.M 2 

K '* l + f e \z\ 2 R r 



(27) 
(28) 



Here, R' is the ratio of the R^/^ values, where R K /^ 
is taken from a single fit to the D° — > K®ir° final state 
only (as opposed to using all three CP— final states un- 

der study), R' = Rk/-k/ R k% ■ and ^ = e si g /ebk g is the 
ratio of the efficiencies of the selection criterion on the 
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helicity angles: f^ K ^ = 0.71 and / e ^o = 0.64. To ap- 
ply these corrections, we first perform a fit of the K®tt° 
final state alone to obtain R K % ■ We then perform the 
simultaneous fit of the CP+ final states, from which we 
take the value of R^/^- Finally, we include the correc- 
tion factors into the CP— final PDF, which will allow the 
likelihood fitter to correctly estimate their influence. The 
parameter \z\ 2 in Eqns. 27 and 28 is extracted from fits 
of the helicity angle distributions in the D° — > K®uj and 
D° -> K° s (p subsamples to the function \z\ 2 + 3 cos 2 9 [29]. 
We subtract the background expected from the Monte 
Carlo simulation, which has been rescaled to match the 



data. We find 



0.065 ± 0.033 in the case of K° 



w, 



and \z\ 2 = 0.217 ± 0.063 in the case of K° s <f>. The un- 
certainties contain propagated uncertainties due to the 
background subtraction. The resulting corrections are: 



,4 true 
A CP(K° s ui) 

,4 true 
A CP(K%4>) 

ptrue 
n K/n{K%u) 

ptrue 



A 



C)1"JH 



cp(k«u)X (1-105 ±0.056), (29) 

x (1.35 ±0.12), (30) 

x (0.9929 ±0.0066), (31) 



/lobs 
^-CPiKlrt,) 



pobs 



= R 



°kUk%*) x (°- 981 ± °' 016 ) ■ ( 32 ) 



In order to assign systematic uncertainties, we propagate 
the uncertainties on the correction factors into the final 
result. 

When calculating Rep through Eq. 5 one has to take 
into account that this equation is an approximation. We 
define the double ratios used to approximate Rcp± as 
R±. They are given by 



R± 



T{B~ ->■ D C p±K-) + T(B+ ->■ D C p±K^ 



T(B- 
T(B- 



D f K-) + T{B+ -> D f K+ 
DfTT-)+T(B+ -^D f n+) 



T{B- -> D C p±tt~) + T(B+ -> D C p±ir+] 



,(33) 



where Df denotes the K ir + final state. These can be 
written as 



R± = 



1 + r 2 B ± 2rs cos 8 b cos 7 

1 + r B r D + 2r B r O COs(SbSd) COS 7 

1 + r 2 B7r r 2 D + 2r B - n r D cos(8 B ^-8 D ) cos 7 
1 + r 2 B7T ± 2r B ir cos 8 b -k cos 7 



,(34) 



where tb-k an d 8b-k arc defined, in analogy to tb and 8b , 



*(<Sb7t-7) - 



A(B~ -> D°n-)/A(B- -)■ D°7r- 



as rg^e 

while r\D and fe are defined as rjje %bD = A(_D° — > 
A"-7r+)/A(L>° -> 7\-7r+). We write Eq. 34 in the 
form R± = Rcp± x (1 + R c ), and we assign a relative 
systematic uncertainty based on the value of the cor- 
rection R c . Taking sin6» c = 0.2257 ± 0.0010 (where 
9c is the Cabibbo angle) and r B =0.104±°;°^ from [3], 
and expressing r D = \V cd V us \/\V ud V cs \ = tan 2 6» c , and 
tbtx = r B ta,n 2 9 c , we find R c w Ar B tan 2 6 C ~ 2.2%. 
Here, we have conservatively assumed values for the co- 
sine terms which maximize R c . We thus assign a relative 



uncertainty of 2.2% to the values of Rcp-, fully correlated 
between Rcp+ and Rcp-- 

We also consider the influence on the measured value of 
Acp of misrcconstructcd signal B candidates, i.e. candi- 
dates reconstructed, in events containing a true B — > DK 
decay with D decaying to the same final state X as the 
reconstructed candidate, from random combinations of 
particles produced in the true B — > DK decay and the 
particles of the ROE. The fraction of these candidates 
ranges from 0.3% to 12% in simulated B ± — s- DcpK^ 
events, depending on the channel. Since we treat this 
component as signal, we implicitly assume that its charge 
asymmetry is equal to the asymmetry in the signal com- 
ponent. We use simulated signal events to estimate the 
ratio between misreconstructed and true B + — i DK + 
candidates and the ratio between misreconstructed and 
true B~ — > DK~ candidates, and find these two quanti- 
ties to differ by less than 0.1%, from which we derive an 
upper limit on the difference between the observed and 
the true value of Acp ■ 

The yield double ratios Rcp± should be corrected by 
the corresponding double ratio of selection efficiencies. 
We find from simulated events that the efficiency double 
ratios are compatible with each other, and their average 
value is very close to unity, (99.46 ±0.23)%. Thus we do 
not correct the central values but conservatively assign 
a relative uncertainty equal to 1 — (0.9946 — 0.0023) = 
0.0077. 

The final PDF doesn't contain an explicit description 
of the conditional parameter A_B sh if t , assuming implic- 
itly that the distribution of A_E shi f t observed in data is 
the same for all the components of the fit. However, the 
distributions are found to be slightly different across the 
components, thus introducing a possible bias in the fit re- 
sults. To estimate the size of this bias, we use simulated 
events to obtain parameterizations of the AE sbitt distri- 
butions of all the fit components and repeat the fits to 
data. We assign the differences compared to the results 
of the nominal fits as the systematic uncertainty. We ex- 
pect this effect to be highly correlated between Acp pa- 
rameters, because the PDFs are similar in each D decay 
channel. Thus they are affected by non-uniform AE s ^±f t 
distributions in a similar way. The same argument holds 
for the Rk/tt parameters. We studied the effect of assign- 
ing a 0%, 50%, and 100% correlation. The uncorrelated 
case gave the largest deviations from the nominal results, 
the fully correlated case gave the smallest. However, the 
variation was found to be at the 10% level. We assign the 
systematic uncertainty corresponding to a correlation of 
50%. 

Table V lists the contributions of the effects discussed 
above. Compared to our previous analysis [9], the sys- 
tematic uncertainty on Aqp+ is reduced due to better 
understanding of the detector intrinsic charge asymme- 
try (the determination of which benefits from the larger 
dataset) and due to improved evaluation of the correla- 
tions among the different sources of systematic uncertain- 
ties. The uncertainty on Acp- is only slightly reduced. 
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By contrast, the systematic uncertainties on Rcp± are in- 
creased due to two additional sources of uncertainty that 
were not considered previously: the bias correction and 
the differences of the AE sii if t distributions among the 
fit components. The systematic correlations between the 
GLW parameters y = (Acp+ , Acp- , Rcp+ , Rcp-) T are 



/ 



C(syst) W\ = 



V 




(35) 



TABLE V: Summary of systematic uncertainties. 



Source 


Acp+ 


A C p- 


Rcp+ 


Rcp- 


Fixed fit parameters 


0.004 


0.005 


0.026 


0.022 


Peaking background 


0.014 


0.005 


0.017 


0.013 


Bias correction 


0.004 


0.004 


0.006 


0.005 


Detector charge asym. 


0.014 


0.014 


- 


- 


Opposite- CP background 


- 


0.003 


- 


0.006 


Rcp± vs. R± 


- 


- 


0.026 


0.023 


Signal self cross-feed 


0.000 


0.001 


- 


- 


s(ir)/e(K) 


- 


- 


0.009 


0.008 


AE shitt PDFs 


0.007 


0.011 


0.029 


0.024 


Total 


0.022 


0.020 


0.051 


0.043 



The statistical correlations among these four quantities 
are: 



1 



C\ 



(stat) [ 



0.0 


-0.08 0.0 


1 


0.0 0.03 




1 0.10 




1 



(40) 



The results are in good agreement with those from our 
previous analysis [9] and the current world averages [21]. 
Figure 2 shows the AE projections of the final fits to 
the CP subsamplcs and Figures 3-5 show tues and T 
projections as well as projections of the fit to the D° —> 
K~tt + flavor mode. 

The statistical significance of a non-zero Acp+ value 
is determined from the maximum value of the likelihood 
function of the nominal fit and that of a dedicated null- 
hypothesis fit, where Aqp+ was fixed to zero, 



•S'stat — V21n(£ nom /£ nu n) — 3.7. 



(41) 



Taking into account systematic uncertainties, the statis- 
tical significance of Acp+ is slightly decreased to: 



■S' 



^81 



stat+syst 



= 3.6. 



(42) 



-y-t 



VIII. RESULTS 

The signal yields returned from the fit for each of the 
D decay mode under study are listed in Table VI. We re- 
construct almost 1000 S ± — >• DcpK^ 1 decays and about 
four times more B^ — > DK^, D — > Kn decays. 



This constitutes evidence for direct CP violation in 
charged B decays and the first evidence of direct CP 
violation in B — > DK . 

We constrain the CKM angle 7, the strong phase 
Sp-, and the amplitude ratio rp from the present mea- 
surement by adopting the frcquentist procedure also ex- 
ploited in [15]. We define a multivariate Gaussian likeli- 
hood function 



TABLE VI: Measured signal yields calculated from the fit re 
suits given in Table VIII using N(B -> DK) = Nf^R 

N(B -)■ Dty) = 
small correlations 



sigO) 



Kot 



K/tt, 

and error propagation neglecting 



D° mode 


N(B ± -> DK^ N{B ± ->■ D^) 


K+K~ 


367 ± 27 4091 ± 70 


tt n~ 


110 ±9 1230 ±41 


K s ir 


338 ± 24 4182 ± 73 


K° s u 


116 ±9 1440 ±45 


kU 


52 ± 4 648 ± 27 


K~n + 


3361 ± 82 44631 ± 232 



£(7,<W\b) 



1 

iV 



cxp 



1 



(y-ytyv-J(y-yt)) (43) 



relating the experimentally measured observables y 
and their statistical and systematic covariance matri- 



ces V c , 



r fit 



I4 ys t with the corresponding truth 



parameters y[ = j/t(7, SB,rp) calculated using Eqns. 3 
and 4. The matrices I4tat and V sys t are constructed from 
Eqns. 35-40. The normalization is N 
We then define a x 2 -function as 



(27r) 2 ^V4; 



X 2 h,S B ,r B ) = -21n£(7,<5 B ,?\B). 



(44) 



The final values of the GLW parameters that we mea- 
sure are: 

A CP+ = 0.25±0.06(stat)±0.02(syst), (36) 

A CP _ = -0.09±0.07(stat)±0.02(syst), (37) 

R C p+ = 1.18±0.09(stat)±0.05(syst), (38) 

R C p- = 1.07±0.08(stat)±0.04(syst). (39) 



Due to the inherent eight-fold ambiguity of the GLW 
method there are eight equivalent minima of the \ 2 ~ 
function, Xmim which correspond to the same value of 
rp and to eight alternative solutions for (7, bp). To eval- 
uate the confidence level of a certain truth parameter (for 
example 7) at a certain value (70) we consider the value 
of the x 2 -function at the new minimum, Xmin(70i ^i r s)i 
satisfying A X 2 = Xm in (7o, S' B , r' B ) - X 2 min > 0. In a purely 



Gaussian situation for the truth parameters the CL is 
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FIG. 2: AS projections of the fits to the data, split into 
subsets of definite CP of the D candidate and charge of 
the B candidate: a) B'^Dcp+K', b) B + ^D C p+K + , c) 
B~^Dcp-K~, d) B + ^D C p-K + . The curves are the full 
PDF (solid, blue), and B^-Dtt (dash-dotted, green) stacked 
on the remaining backgrounds (dotted, purple). The region 
between the solid and the dash-dotted lines represents the 
B^rDK contribution. We show the subsets of the data sam- 
ple in which the track h from the B decay is identified as a 
kaon. We require candidates to lie inside the signal-enriched 
region defined in Sec. IV, except for the plotted variable. 



FIG. 3: jties projections of the fits to the data, split into 
subsets of definite CP of the D candidate and charge of the 
B candidate: a) B~ -»■ D C p+K~ , b) B + -* D C p+K + , c) 
B~ -} Dcp-K~ , d) B + — >• Dcp~K + . We show the subsets 
of the data sample in which the track h from the B decay is 
identified as a kaon. See caption of Fig. 2 for line definitions. 
Only a subrange of the whole fit range is shown in order to 
provide a closer view of the signal peak. 



given by the probability that A% 2 is exceeded for a x 2 
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FIG. 4: T projections of the fits to the data, split into sub- 
sets of definite CP of the D candidate and charge of the 
B candidate: a) B~ -* D C p+K' , b) B + -> D C p+K+, c) 
B~ ->■ D C p-K~ , d) B + -> Dcp-^ + . We show the subsets 
of the data sample in which the track h from the B decay is 
identified as a kaon. See caption of Fig. 2 for line definitions. 



distribution with one degree of freedom: 



1-CL = 



V2r(i/2) 



e -t/2 t -l/2 dt (45) 
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FIG. 5: Projections of (a) AE, (b) tties, and (c) J- variables 
of the fit to the B ± -> D/i ± , D° -> ^"tt" 1 " flavor mode. No 
requirements are put on the PID of the track h from the -B 
decay and on the fit variables not plotted. See caption of 
Fig. 2 for line definitions. 



A more accurate approach is to take into account the 
non-linearity of the GLW relations, Eqns. 3 and 4. In 
this case one should consider A^ 2 as a test statistic, and 
calculate (1 — CL) by means of a Monte Carlo procedure, 
described in the following. For a certain value of interest 
(7o), we: 



1. calculate A% 2 



K 2 ■ 
Amin 



ho,8' B ,', 



x 2 - 

A mm 



as before; 



Ax 2 



2. generate a "toy" result A' CP ±, R' C p±, using Eq. 43 
with values 70, S' B , r' B as the PDF; 

3. calculate A% 2 ' of the toy result as in the first step, 
i. e. minimize again with respect to 5 B and r B ; 
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4. calculate (1 — CL) as the fraction of toy results 
which perform better than the measured data, i.e. 
1 - CL = N(A X 2 > A x 2 ')/N t oy. 

Figures 6 and 7 illustrate 1-CL as a function of 7 and tb 
as obtained from this study. From these distributions we 
extract 68% and 95% CL confidence intervals for 7 and 
tb , as summarized in Table VII. Due to the 7 «-> 8b am- 
biguity of the GLW method, the ID CL intervals for 5b 
are identical to those for 7. At the 68% CL we are able to 
distinguish six out of eight solutions for 7 (and 5b), two 
of which are in good agreement with the current world 
averages [25]. At the 95% CL we are able to exclude the 
intervals [0°,7.0°], [173.0°, 187.0°] and [353.0°, 360°] for 
7 and 5b- For tb we deduce at 68% CL: 



tb 



0.35to;li(stat + syst). 



(46) 



TABLE VII: 68% and 95% CL intervals for the parameters 7, 
and tb, taking into account both statistical and systematic 
uncertainties. The confidence intervals for 5b are identical to 
those for 7 due to the intrinsic 7 •<-> 5b ambiguity of the GLW 
method. 



7 



mod 180 I 



Tb 



68% CL 



95% CL 



[11.3,22.7] 

[80.8, 99.2] 

[157.3, 168.7] 

[7.0, 173.0] 



[0.24, 0.45] 



[0.06,0.51] 



In order to facilitate the future combination of these 
measurements with the results of the Dalitz plot analysis 
of B ± -> DA'i, D -► K°h+h- decays (h = n,K) [16], 
we recompute the GLW parameters after excluding from 
the nominal fit the D C p- -> K°(f> ((f) — > K + K~) sub- 
sample. The sample obtained in this way is statistically 
independent of that selected in [16]. The final values of 
the GLW parameters that we measure in this case are: 

A CP+ = 0.25±0.06(stat)±0.02(syst), (47) 

A CP - = -0.08±0.07(stat)±0.02(syst), (48) 

R C p+ = 1.18±0.09(stat)±0.05(syst), (49) 

R C p- = 1.03±0.09(stat)±0.04(syst). (50) 

The statistical correlations among these four quantities 
arc: 



a 



(stat) 



[y\ = 



1 -0.08 
1 0.04 
1 0.09 



(51) 



and the systematic correlations are 



a 



1 0.56 -0.06 
1 

(syst)[J/J = 1 0.12 

1 



(52) 
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FIG. 6: 1-CL as a function of 7 (top) and rs (bottom). 
Both statistical and systematic uncertainties are taken into 
account. For the angle 7, the plot is identical in the range 
[180°, 360°]. The horizontal lines show the 68% CL (dashed) 
and the 95% CL (dotted). Due to the symmetry of Eqns. 3 
and 4 the plot for the strong phase 5b is identical to the one 
for 7. 



To compare the results obtained after removing the 
Dcp- — > K®<f> subsample with those from the B ± — > 
DK^jD — > K®h + h~ analyses, which are expressed in 
terms of the variables x± = tb cos(<5b ± 7) and y± = 
r_Bsin(i5_B ±7), we use the GLW parameters measured 
in this way to determine the quantities x± through the 
relations: 



%± = J [Rcp+i^-T A CP+ 
We obtain 



R CP -{lT A C P-)\ . (53) 



x + = -0.057 ±0.039(stat)±0.015(syst), (54) 
x_ = 0.132 ±0.042(stat)±0.018(syst). (55) 

These results are in good agreement with the current 
world averages [21] and have precision close to the single 
most precise measurements [16]. We also measure rfj, 
which provides a constraint on x± and y± via r B = x± + 
j/|, from 



1 



[Rcp+ + Rcp- — 2) . 



(56) 



We determine: 



0.105 ± 0.067(stat) ± 0.035(syst) . (57) 
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FIG. 7: Contours at 68% (dotted, red) and 95% (solid, green) 
2-dimensional CL in the (7, rs) and (7,5s) planes. See also 
the caption of Fig. 6 regarding symmetries. 



The constraints that could be placed on the quantities 
y± from these measurements, by exploiting the relation 
r B = x\ ± y± , are much weaker than those provided by 
the B ± -> DK±, D -> K°h + h~ analysis. 

As a final check of consistency we consider the quantity 



— Acp+Rcp+ + Acp-Rcp- 



(58) 



From Eqns. 3 and 4 one expects a to satisfy a = 0. We 
measure a = 0.19±0.11(stat + syst), which is compatible 
with 0. 



IX. SUMMARY 

Using the entire dataset collected by BABAR at the 
e + e~ center-of-mass energy close to the T(4S) mass, we 
have reconstructed B^ — > DK ± decays, with D mesons 
decaying to non-CP (Ktv), CP-even (K + K~ , ir + ir~) and 
CP-odd (K°n°, K% K°lo) eigenstates. 

Through an improved analysis method compared to 
the previous BABAR measurement [9] and through an en- 
larged dataset, corresponding to an increase in integrated 



luminosity at the T{AS) peak from 3481b" 1 to 426 fb _1 , 
we obtain the most precise measurements of the GLW 
parameters Acp± and Rcp± to date: 

A CP + = 0.25±0.06(stat)±0.02(syst), 

A C p- = -0.09±0.07(stat)±0.02(syst), 

R C p+ = 1.18±0.09(stat)±0.05(syst), 

R C p- = 1.07±0.08(stat)±0.04(syst). 

We measure a value of Aqp+ which is 3.6 standard 
deviations from zero, which constitutes the first evidence 
for direct CP violation in B — > DK decays. 

From the measured values of the GLW parameters, we 
extract confidence intervals for the CKM angle 7, the 
strong phase Sp , and the amplitude ratio tb , using a fre- 
qucntist approach, taking into account both statistical 
and systematic uncertainties. At the 68% CL we find that 
both 7 and dp (modulo 180°) belong to one of the three 
intervals [11.3°, 22.7°], [80.8°, 99.2°] or [157.3°, 168.7°], 
and that 

r B e [0.24,0.45]. 

At 95% CL, we exclude the intervals [0°,7.0°], 
[173.0°, 187.0°] and [353.0°, 360°] for 7 and S B , and mea- 
sure 

r B e [0.06,0.51]. 

Our results are in agreement with the current world av- 
erages [25]. 

To facilitate the combination of these measurements 
with the results of our Dalitz plot analysis of B ± — s- 
DK±, D -> K° s h + h- (h = K, n) [16], we exclude the 
D —> K®4>, 4> — >• K + K~ channel from this analysis - thus 
removing events selected also in [16] - and then deter- 
mine 

A C p- = -0.08±0.07(stat)±0.02(syst), 
R C p- = 1.03±0.09(stat)±0.04(syst). 

For comparison with the results of the B ± — > DK^ , D — > 
K®h + h~ analyses, which are expressed in terms of the 
variables x± = r B cos(S B ± 7) and y± = r B sm(S B ± 7), 
we express our results for the GLW observables in terms 
of x+ and X-. We measure 

x+ = -0.057 ±0.039(stat)±0.015(syst), 
x_ = 0.132 ±0.042(stat)±0.018(syst), 

at 68% CL. These results arc in good agreement with the 
current world averages [21] and have precision compara- 
ble to the single most precise measurements [16]. We also 
evaluate r B after the exclusion of the D — > K®<fi channel, 
and obtain a weak constraint on r B = x± ± y±: 



0.105 ± 0.067(stat) ± 0.035(syst) 



at 68% CL. 
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TABLE VIII: Fit result of the three final fits to data, before correcting for fit biases (see 


Section VII). 


Parameter 


KK 


7T7T 


K S TT 


K%u 


kU 


K-k 


4 sig(K) 
A CP 


0.242 ± 0.065 




-0.089 ± 0.066 




-0.008 ± 0.022 


/isigU) 
"-CP 


0.003 ± 0.015 




-0.009 ± 0.014 




-0.0116 ±0.0050 


Rk/tt 


0.0897 ±0.0063 




0.0808 ± 0.0056 




0.0753 ± 0.0018 


m 


0.0204 ± 0.0030 




0.0206 ± 0.0029 




0.02143 ±0.00089 


A BB 
A CPJ 


-0.004 ± 0.045 


-0.043 ± 0.047 


n/a 


n/a 


n/a 


-0.043 ±0.017 


A qq 
CP,p 


0.012 ±0.016 


-0.016 ±0.018 


-0.002 ±0.011 


0.012 ±0.020 


-0.069 ± 0.060 


-0.027 ±0.016 


All 
CPJ 


-0.004 ±0.011 


-0.0044 ± 0.0098 


0.0021 ± 0.0071 


-0.004 ±0.013 


0.001 ±0.039 


-0.0016 ±0.0068 


ff 


0.326 ± 0.026 


0.49 (fixed) 


0.520 ± 0.030 


0.27 (fixed) 


n/a 


0.396 ±0.018 


°%,1 


0.160 ±0.016 


0.258 ± 0.023 


0.206 ± 0.014 


0.175 ±0.034 


0.2758 ± 0.0092 


0.198 ±0.014 


CT qq 


0.1742 ±0.0020 


0.2047 ± 0.0024 


0.1546 ±0.0015 


0.1963 ±0.0028 


n/a 


0.1965 ±0.0017 


rr qq 


0.312 ±0.011 


0.329 ±0.011 


0.3541 ± 0.0068 


0.317 ±0.019 


0.447 ±0.014 


0.3068 ± 0.0061 


° T,r,1 


0.231 ±0.014 


0.268 ±0.018 


0.275 ± 0.020 


0.238 ±0.013 


n/a 


0.237 ±0.010 


„BB 

a AE,mm 


n/a 


n/a 


n/a 


n/a 


n/a 


0.01048 ± 0.00057 


n qq 

AE 


-0.96 ±0.14 


-0.71 ±0.14 


-0.924 ± 0.099 


-1.04 ±0.18 


-0.48 (fixed) 


-0.88 ±0.10 


VAE 


-2.62 ±0.32 


-1.36 ±0.57 


-1.80 ±0.35 


-2.87 ±0.59 


-0.95 ±0.75 


-1.527 ±0.092 


a AE 


16.63 ±0.27 


14.82 ± 0.49 


17.01 ±0.29 


16.10 ±0.52 


15.82 ±0.60 


15.424 ± 0.076 


MmES 


5278.56 ± 0.12 


5278.61 ± 0.20 


5278.62 ±0.12 


5278.50 ± 0.22 


5278.99 ± 0.25 


5278.586 ± 0.033 


"•jTiES,! 


2.207 ±0.081 


2.12 ±0.15 


2.299 ± 0.084 


2.11 ±0.16 


2.33 ±0.17 


2.210 ±0.022 


&rriES,r 


2.897 ±0.081 


2.83 ±0.15 


2.922 ± 0.084 


3.08 ±0.16 


2.72 ±0.17 


2.852 ±0.023 


N* B 


79 ±29 


346 ± 52 


176 ± 43 


180 ± 48 


3 (fixed) 


328 ± 40 


N fB 


1430 ± 82 


1517 ±142 


1930 ± 102 


1195 ± 109 


119 ±20 


7717 ± 170 


N? 


4005 ± 69 


3456 ± 76 


8587 ± 101 


2675 ± 68 


284 ± 17 


4722 ± 77 


Nf 


10890 ± 125 


13019 ± 176 


21657 ± 172 


6673 ± 124 


716 ± 29 


28007 ± 205 


<? w 


4091 ± 70 


1230 ±41 


4182 ± 73 


1440 ± 45 


648 ± 27 


44631 ± 232 



